A Simple Predator-Prey Population Model with Rich Dynamics

A non-smooth switched harvest on predators is introduced into a simple predator-prey model with logistical growth of the prey and a bilinear functional response. If the density of the predator is below a switched value, the harvesting rate is linear; otherwise, it is constant. The model links the well studied predator-prey model with constant harvesting to that with a proportional harvesting rate. It is shown that when the net reproductive number for the predator is greater than unity, the system is permanent and there may exist multiple positive equilibria due to the effects of the switched harvest, a saddle-node bifurcation, a limit cycle, and the coexistence of a stable equilibrium and a unstable circled inside limit cycle and a stable circled outside limit cycle. When the net reproductive number is less than unity, a backward bifurcation from a positive equilibrium occurs, which implies that the stable predator-extinct equilibrium may coexist with two coexistence equilibria. In this situation, reducing the net reproductive number to less than unity is not enough to enable the predator to go extinct. Numerical simulations are provided to illustrate the theoretical results. It seems that the model possesses new complex dynamics compared to the existing harvesting models.